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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE NO. 186^ 


EXTENDED APPLICATIONS OF THE HOT-WIRE ANEMOMETER 
By Stanley Corrsin 


SUMMARY 


Two new fields of application of the hot-wire anemometer are 
proposed, and the appropriate response equations and measuring procedures 
are developed. 

The first analysis leads to a method for the measurement of 
physically significant statistical quantities in a turbulent flow with 
heat transfer j for example, the turbulence levels, the temperature 
fluctuation level, the turbulent heat-transfer coefficient, the velocity 
scale, the temperature scale, and some spectrums- 

The second analysis involves the use of the hot-wire in the turbulent 
isothermal mixing of two appropriately different gases- If the thermal 
conductivity of the mixture is known and is a monotonic function of the 
relative concentration, it is possible to measure the mean velocity and 
mean concentration at any point. If no data are available on the thermal 
conductivity of the mixture, this additional unknown can be determined by 
an additional measurement. Furthermore, it is also possible to measure 
the various statistical functions of the fluctuating velocities and the 
local concentration fluctuation, provided, again, that the ther mal 
conductivity is a known monotonic function of the concentration. 

Although the details of the present analysis are dependent upon 
the accuracy of King's equation for the rate of heat loss from fine wires, 
the general approach is equally valid for any possibly more accurate 
equation that may be deduced. 

\ 


INTRODUCTION 


The hot-wire anemometer has been used for many years in the 
measurement of the mean velocities of flowing gases (reference l) and, 
more recently, in the measurement of various statistical q uan tities 
associated with the random velocity fluctuations occurring in turbulent 
flow (references 2 to 7 and others) • It has also been used as a micro- 
phone (reference 8) . This anemometer consists essentially of an 
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electrically heated wire of extremely small mass, the temperature of 
which varies in response to rapid changes in the instantaneous flow 
(cooling) velocity. If the heating current is maintained constant, the 
resulting voltage fluctuations, the intensity of which is expressible 
as a function of the velocity fluctuations, may he amplified and measured. 
Attenuation and phase lag, due to the nonzero heat capacity of the wire 
and the finite rate of electrical heating, may he compensated hy a 
suitable network between two of the amplifier stages. The essentials 
of the technique are described in reference 2. 

In a detailed consideration of turbulent fluid flow, there are 
several statistical functions of the randomly fluctuating velocities 
that are of importance. The simplest of these are the following: 

u' 

(a) The three components of the turbulence level or intensity — , 

v ' w 1 ^ 

and =r which, squared and added, give the ratio of kinetic energy of 

random fluctuation to directed kinetic energy. 

(b) The double correlation functions, that is, the statistical 
correlation between two components at separate points; for 

u 1 u 2 

example, R P = • Von Karman (reference 9) and Von Kar man and 

u ' u ' 

1 2 

Howarth (reference 10) have shewn that only two of these functions are 
independent for isotropic turbulence. 

(c) The energy spectrum of the turbulent fluctuations. Taylor 
(reference 11) has demonstrated that this is the Fourier transform of 
one of the two principal components of the correlation tensor.^ - 


(d) The scale of turbulence, defined by Taylor (reference 12) for 
isotropic turbulence as the integral of one of the double correlation 

^ 00 

functions; for example, L = Rg^dy* 

J 0 


(e) The microscale of turbulence, identified with the rate of 
dissipation of turbulent energy into heat, the size of the smallest 
eddies, was also defined by Taylor (reference 12), who showed this to 
be the abscissa intercept of the vertex-tangent parabola of the R 2 (y) 

curve in isotropic turbulence. 


1 


In low-intensity isotropic turbulence. 
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(f) The turbulent shear stress t = -puv was first identified by 
Osborne Reynolds in his classic formulation of the equations of motion 
for turbulent flow (reference 13)* This stress is, of course, a measure 
of the lateral rate of momentum transfer or diffusion. 

All the foregoing functions have been measured for various turbulent 
flows (references 6, 1 , 14, 15, 1 6, 11 , 18 , 19, and others), principally 
with the hot-wire technique, although the most detailed measurements have 
of necessity been confined to investigations of the "isotropic turbulence 
far behind- grids placed in a uniform air stream. 

There are other kinematic quantities of considerable interest: 

(g) Batchelor and Townsend (reference 20) have recently measured 
three statistical functions of the time derivatives of the velocity 

, and ( . These enter into terms 

\dt/ 

describing the time rate of change of vorticity in isotropic turbulence. 

(h) The rate of lateral kinematic diffusion of turbulent energy 

in a turbulent shear flow (~ pu^v) has apparently not been measured 
yet, although suitable techniques are fairly well known. 

(i) The pressure fluctuation diffusion term (~ pv) in the 
turbulent energy equation has not been measured, principally because of 
the difficulty of 'obtaining an instrument that responds to static-pressure 
fluctuations independently of velocity fluctuations. 

In a detailed consideration of turbulent motion, it is of interest 
to obtain information on the heat and material diffusion coefficients as 
well as the momentum and energy diffusion already mentioned. In fact, 
this information should be useful both as an end in itself and as 
additional evidence to be applied in obtaining a better understanding of 
the nature of turbulence. Consequently, it seemed worthwhile to investi- 
gate the possibilities of exploiting the hot-wire beyond the simple 
anemometric applications. 

The rate of cooling of a body immersed in a fluid stream clearly 
depends upon: (a) The stream velocity, (b) the temperature difference, 

(c) the physical constants of the fluid, and (d) the physical constants 
of the wire. For the simple hot-wire anemometer, conditions (b), (c), 
and (d) are effectively fixed, and the instrument response measures 
condition (a). Changes in stream temperature certainly involve changes 
in condition (b) and possibly in condition (c). Thus the instrument 
responds to temperature changes or fluctuations. Changes in the 


fluctuation 
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relative concentration of a flowing mixture of two fluids with different 
physical constants lead to variations in condition (c). In principle, 
the hot-wire instrument responds to concentration changes or fluctuations. 

Considerations of algebraic and experimental complication have 
restricted the present treatment to cases in which fluctuations in 
velocity and only one of the aforementioned two quantities coexist, 
although this limitation is not dictated by principle. 

With this experimental technique, the new quantities which’ it may 
be of interest to measure are the temperature fluctuation level 

the spectrum of temperature fluctuation, the heat diffusion (~ $v ) , the 
correlation coefficient I 1, and the thermal scale and microscale, 

\e 2 J 

analogous to Taylor's kinematic scales, as well as the corresponding 
statistical functions of the concentration fluctuation in the turbulent 
mixing of two gases. 

The fact that direct measurement of dv should be possible with a 
hot-wire instrument has been known for many years (see, for example, 
remarks in reference 2l) , but no actual work appears to have been 
carried out. 

Distributions *of momentum, heat, and material transport coefficients 
in a shear flow are directly calculable from measurements of the mean 
velocity, temperature, and concentration profiles, respectively. However, 
it has been found that the mean profile is relatively insensitive to 
phenomenological assumptions on the transfer coefficient (references 16 
and 19 ), so that direct experimental determination of the transfer is 
desirable. 



SYMBOLS 


U Instantaneous total velocity 

U mean velocity 

u instantaneous velocity fluctuation in direction of 

mean velocity (U - U) 


normal components of velocity fluctuation 
root -mean- square values of u, v, and v 
instantaneous absolute temperature of gas 

mean absolute temperature 

instantaneous temperature fluctuation ^T a - T a 
root -mean- square value of -Q 

absolute temperature of reference medium (e • g . , 
the air at rest, for a free Jet) 


instantaneous concentration of secondary gas j in 
present applications, 

_ Hols per unit volume of the secondary gas 
Total mols, of air and gas, per unit volume 

mean T 

instantaneous concentration fluctuation (T - f ) 
root -mean-square value of 7 

instantaneous resistance of electrically heated 
hot-wire (unless otherwise noted) 

mean R 


instantaneous resistance of unheated hot-wire 
mean R a /resistance at 
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instantaneous "equilibrium" hot-wire resistance; 
that is, resistance hot-wire would have at any 
moment if it. had zero lag 

mean B e 

resistance of wire at T r 
resistance of wire at 0° C 
heating current 

voltage fluctuation across hot-wire 
hot-wire length 
hot-wire diameter 

specific resistivity of wire material 


thermal coefficient of change of resistance of 
wire material 

v sensitivity of directional meter 

coefficients in assumed parabola relating thermal 
conductivity to concentration in a gas 
mixture 

7 sensitivity of a simple hot-wire 
ohms (electrical resistance) 
density of wire material 


mass of wire 

specific heat of wire material 
thermal conductivity of fluid at T a 
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*8 

^ = k g " k a 


C Pa 

c % 

Ao p ■ °p g ■ °p 0 


thermal conductivity of fluid at T a 
thermal conductivity of air 
thermal conductivity of secondary gas 

constant-pressure specific heat of fluid at T a 
constant -pressure specific heat of fluid at To 

a. 

constant -pressure specific heat of air 
constant -pres sure specific heat of secondary gas 


P 


P 


p a 

p g 

k', c p ', p' 

= Pg ' Pa 

H 


T 

t 


density of fluid at T a 
density of fluid at "T a 
density of air 
density of secondary gas 
fluctuations 

heat energy in hot-wire 
wire temperature 
time 
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CONVENTIONAL HOT -WIRE -ANEMOMETER RESPONSE EQUATIONS 


Following the approach of Dry den and Kuethe (reference 2) , King's 
semi empirical form is assumed for the rate of heat loss from a wire 
immersed in a flowing fluid. For the time rate of increase of heat 
energy in the wire this gives 

“ = i 2 R - (A' + B' n/u)(t - Tg) (1) 

where the physical make-up of A' and B' as deduced hy King is 

A' = alk 

B 1 = hi ^dcppk 


Here a and b are regarded as numerical constants, although they are 
dimensional, and should show some variation if a wide enough temperature 
range is considered. 

The term A'(T - T Q ) presumably represents heat loss due to free 

convection and radiation. Therefore, it is clear that the constant a 
must include such physical quantities as the acceleration of gravity, 
the thermal expansion coefficient of the fluid, the specific heat of the 
fluid, a temperature function, a radiation constant, and so forth. The 

term B 1 ^U(T - T a ) represents forced convection, and, after k has been 

factored out, the term can be written as proportional to the product of 
a Reynolds number and a Prandtl number. 

Equation (l) can be written as 


4.2ms dR 

R Q a dt 



(la) 


and if equilibrium conditions are considered such that — =0, writing 

dt 

the equilibrium value of R as R Q , the following equation is obtained: 


i He 


®e " ®a 


= A + B\[u 


(2) 
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where 



In considerations of averaged or steady-state operation, R 0 in 
equation (2) may he replaced by R, and this equation is the usual mean 
velocity calibration of the hot-wire. 

King deduced his equation for steady-state operation, and a 
primary assumption of the hot -wire -anemometer theory is that this rate 
of heat loss is independent of acceleration. 

A convenient alternative form of equation (2) is 

- 1 - — = A - i 2 + B \|u (2a) 

R - R n 


-• For the measurement of turbulence, interest is centered in the 
voltage fluctuation set up across the hot-wire due to a small fluctuation 
in velocity. The analysis for large fluctuations has proved intractable, 
but, as illustrated in appendix A of reference 22, the error committed in 
applying the small perturbation results in the measurement of rather large 
fluctuation levels is not necessarily excessive. 

Letting 

U = U + u 
R = R + r 
e = ir 

where 

•=■ « 1 and r 

u ■ s - s a 


« 1 
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and substituting these into equation (2) 


yields 


i 2 (B + r) 




Keeping only the linear terms of the appropriate binomial expansions, 
otherwise neglecting the squares of small quantities, and substituting 
from equation (2) for (A + B \fu) , there is finally obtained 



[i 2 S - a(r - R a )] £ 


(3) 


Thus, the coefficient of is the sensitivity of a single hot-wire to 

u 

velocity fluctuations along the main stream direction* 

Since the wire has a nonzero heat capacity and only a finite rate of 
electrical heating, there is a time lag in the response which normally 
leads to an appreciable drop in the curve of hot-wire response against 
frequency above some particular frequency. The characteristic "time 
constant M can be computed approximately by substituting from equation (2) 
into equation (la) to obtain 

4.2ms dR * R a (R e ~ ( 4 ) 

E Q a dt R e - R fl 

Dryden and Kuethe (reference 2) transform this further but find no 
useful general solution. Therefore, they use the assumption of smal 1 
perturbations to justify the replacement of (R e - R a ) by (R - R a ) in 

the denominator of the right side, which leads to the simple form 


dR 

dt 


+ 


M 


M 


(5) 


or 


S< K - S > + 5< B - E e> = 0 


where the time constant 


M = 


4.2ms(R - R a ) 


2- 
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is the time required for (R - R) to teccane equal to — times an 

e 

original step-function difference, say, (R^ - R) . The definition of M 
is, of course, not dependent upon the assumption of small perturbation. 
Further details are discussed in reference 2. 


Physically, the time constant can be considered as proportional to 
the product of three factors, functions of the wire material, size, and 
operating conditions, respectively: 



( 6 ) 


The attenuation and phase lag in the hot-wire response are compensated 
up to appropriate frequencies by a suitable network between two of the 
amplifier stages- 

The technique for the measurement of specific statistical kinematic 
quantities in a turbulent flow need not be discussed in detail here- 
Adequate descriptions are given in other papers (references 2, 7, 18, and 

„ i 

others) 

a measurement of 


u 


It is sufficient to point out that =- is computed directly from 




as given in equation (3) 


The components r— > 

U U 

and uv are measured with directionally sensitive hot-wire instruments, 
customarily in the 'form of an X, with the polar axis perpendicular to the 
mean velocity direction. Each of the two wires in such an instrument, 
being inclined to the mean flow, responds to both longitudinal and lateral 
fluctuations, so that, instead of equation (3), there results 


P u , v 
e-, 0 = &=- ± e— 
U U 


for the two wires. Thus, for example. 


ir ~v (e i + 62)2 


F v (e i ■ e 2 )2 


uv 2 p 

— ~ e. - 


w O-j 

5 s 
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The correlation between u fluctuations at two different points is 
obtained in the same way as the shear. Integration of the correlation 
curve gives the scale of turbulence. 

The energy spectrum of turbulence is obtained conventionally with 
a narrow band-pass filter. 

Some corrections for the effect of wire length have been worked out 
in reference 7, but further detailed analysis along that line is necessary, 
particularly for the measurement of correlation functions and shear and 
turbulence spectrums. K. Zebb (reference 32) has made some approximate 
computations for the length correction of X meters in v' measurements. 


MEASUREMENT OF MEAN VELOCITY AND TEMPERATURE 


In applications of King's equation for thermal equilibrium 

i 2 R 

* - »a 


A + B 


{s 


where 



B 


bZ 

R^a 

0 


4 


dCpPk 


to the measurement of mean velocity and temperature in flow with a 
temperature gradient, the_only change from the first part of the preceding 
section is that A and B now vary from point to point in the flow. 

These variations are determined from the changes in the physical constants 
of the fluid. The present investigation is restricted to atmospheric air. 

Since air is very nearly a perfect gas. 


_P 

Pr 



(7) 


where the subscript r corresponds to an appropriate reference temperature. 
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From the International Critical Tables (reference 23), the thermal 
conductivity of air can be approximated for temperatures up to a few 
hundred degrees centigrade by the empirical equation 


k_ _ T r + 12 2 / Ta\3/ 2 
k r “ T a + 125 


This can be approximated by a straight line over a fairly wide temperature 
range, and, in fact, the linear approximation is used in the next section. 

The variation in Cp with temperature can be neglected in the 
present discussion, from 0° C to about 300° C- Thus, 


A_ _ k_ 
^ " kr 



In practice, it is convenient to have these expressions in terms of 
wire resistances. With 


K a = R 0 [l + a(T a - 273)] 
K r = Rq [l + a(T r - 273)] 


( 10 ) 



equation ( 9 ) can toe written in the form 



These are presented in figure 1 for T r = 293° atos. 
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With these values of A and B, the mean velocity is determined_ 
just as in the case of the simple anemometer. The mean temperature T a 

is, of course, obtained directly from equation (10) which can he combined 
to give 


T 


a 


= T r 




(13) 


which, for a platinum wire becomes approximately 

T =t(^] 

a r \R r J 


It should be emphasized that the accuracy of these results, as well 
as all those in the following sections, depends upon the accuracy of the 
physical form of the cooling terms in King's equations. The linear 

i 2 R 

variation of — with 

R - R- 




has been checked with considerable 


accuracy by many experimenters; figure 2 is a typical isothermal mean- 
velocity calibration of a hot-wire anemometer. However, there has 
apparently been no attempt to check the exponents of the other physical 
quantities as derived by King. Appendix A contains a discussion of this 
problem, as well as some preliminary experimental results. 


RESPONSE OF A HOT-WIRE IN AN AIR STREAM WITH 
VELOCITY AND TEMPERATURE FLUCTUATIONS 
Voltage Fluctuation 


In order to deduce an equation analogous to equation ( 3 ) , start 
with King's equation for thermal equilibrium. 




1/2 
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With small fluctuations in velocity and air temperature, there are 
introduced 


R = R + r 


®a ^a + r a 


U = U + u 
P = P + p' 


k = k + k ' 


c p “ q P 


so that King's equation becomes, -writing 


al 

R 0 a 


= P and 


h l ydc 
R 0 a 




i 2 (E + r) 


(R - R a ) + (r - r a ) 


= P(k + k') + Q Qp + p')(k + k')(u + u)J 


1/2 


Using a "binomial expansion for the denominator of the left side 
and neglecting products of small quantities on both sides, 


fix + - 5 > r ) ■ PS + k ’> 


+ Q(pkU + p'kU + pk'U + pku) 


Keeping the linear terms of the binomial expansion of the last part, 


i 2 R 

— + 


E ' R a (R - R a ) 


i 2 JC- _ X 

r-g (Br a ‘ W = p ( k + k ') 

p - p v * 


«(Su) l/2 [r + 0 * % * §) 


In order to get the fluctuation equation, subtract the averaged King's 
equation, 


i 2 R 


R - R 


— = Pk + Q(pkU) 


1/2 


a 
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which leaves 


■ v) = pk ' + 2 - <pku)1/2 (5 + f + 1) ■ (i « 


and the next step is to express r a , k 1 , and p 1 in terms of the sm&l 1 
temperature fluctuation d . 

Obviously, 


r a = R Q ad 


(15) 


As mentioned in the previous section, the rather complicated empirical 
expression for the temperature dependence of the thermal conductivity of 
air (equation (8)) can be well approximated over a large temperature, range 
by a single straight line- Hence, with considerable accuracy, there can 
be written 


k' = nd (16) 

where n is sensibly constant over a range of a few hundred degrees 
centigrade. 


For the density fluctuation, again air is considered as a perfect 
gas and the p(T) hyperbola in the vicinity of T & is approximated by 
its tangent. Thus, 



P + P 1 


Therefore, 


P + P* 


_ Ta 

P = 

T a + fl 



and for 


■=- << 1, there results 
T a 


P' = 


d 



(17) 
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Substituting equations (15), (16), and (17) and introducing the 
hot-wire voltage fluctuation e = ir, as well as the expressions 
for P and Q, 


i 2 KRoa 


1E« 


"0 - 


(E - (R - R a ) 2 


= + it(4c p Sd) 1/ 2 (- i , n i , i 1 ] ( 18 ) 

V V V ' V T * k v) 


Calculation with typical, numerical values (see appendix B) shows that, 
for most practical purposes, the sum of the two d-terms . in the final 
parenthesis of equation (l8) can be neglected in comparison with the first 
term in the equation, so they will be omitted in the rest of the discussion. 

Collecting the remaining terms. 


^a 

i 2 RR~a , 

0 aZn 

d + bI ( 

(R - R a ) 2 

(H - R a ) 2 E 0 a 

2R 0 a\ 

The functional form desired is e = 

e(t, ij. 



\ e V 


R 

a - R r 


9 = T - T .= - 
a x r 




R 0 a 


U 


the d-term is multiplied by ( tr) • It is also convenient to 


R 0 a ej 


replace ^ ^(dCpPkTj) ' by its equivalent as given in King's equation. 


£L/^„ ^i/2 

RqO 

Then, finally. 


e ~ 


E a 

iRo 


T 

— V 4 


2- < E “ E a) 2 

i R - aZn — 

d R - R»r 2 - _ 1 

- " - il-A E- R a 

( E ° a )" 

9 2iR L a 

ci 


u 


(19) 


which obviously reduces to equation ( 3 ) when the temperature fluctuation 
level becomes much smaller than the velocity fluctuation level. The value 
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of A Is determined from the mean-speed calibration of the hot wire, 

, aZ A 

and, of course, = — • 

RqQ, k 


Time Constant 

The general equation is given under CONVENTIONAL HOT -WIRE -ANEMOMETER 
RESPONSE EQUATIONS as follows: 

- ( A + B Vu)( K -E a) (la) 


Equation (2) describes a hypothetical equilibrium reached at the 
instantaneous flow conditions. Substituting for (A + B >/u) from equation (2), 
there is obtained 

4.2ms dR _ i 2 R a (Ne ~ K) (20) 

Ro a R e " R a 


which is the equivalent of equation (4) of reference 2, with the essential 
difference that here R a is a fluctuating function of time- 

Since, as for the simple anemometer, the general solution of equation (20) 
is apparently not expressible in any applicable form, possible simplifica- 
tions are now considered for the restriction of small fluctuations, that is. 



« 1 


and 



« 1 


Hence R a may be replaced by R a in the numerator of the right side 
of equation (20) • The denominator may be examined in somewhat the same 
way as Dry den and Kuethe obtained their equation (15) j that is, there may 
be written 

E e - R a = ( R " *0 + ( R e " R ) + (*a ’ ^ 

Then it is seen that the last two parentheses are small compared with the 
first if appreciable heating current is used. In some cases, when only 
the temperature fluctuation level is being measured, relatively little 
heating current is used. With negligible sensitivity to velocity, it is 
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apparent that R e - R a = Constant and, therefore, may be replaced by the 
average R g - R a . However, for small fluctuations, R e = R, so that 
R e “ R a ma y replaced by R - R a . Equation (20) becomes approximately, 

4.2ms dR = ^a( R e ~ R ) (2l) 

Roa dt R - R a 


or 


, dR 

M S ■ S e 


R 


(21a) 


where 


4.2ms(R - Ra) 
i 2 RaR 0 a 


( 22 ) 


so that equation (21a) is identical with Dryden and Kuethe's equation ( 15 ), 
and M 1 = M, the simple anemometric time constant. Clearly, from the 
foregoing analysis, this will not be true for large fluctuations. 

With negligible heating current, equation (la) becomes 

^ ® - -( A + B \AS)(S - »a) 

RqCC» CLu 

From equilibrium, R e = R a - Therefore, 

(A + B \/u) (R e - r) 

^ = R - R (21a) 

dt 6 

4.2ms 

R Q a(A + B \fd) 


or 


where 


4.2ms dR 
Roa. dt 


M' 


M' = 


(23) 
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But 


(A + B \ju) = 


1 2 R 
R - R 


a 


R - R c 


for the case of extremely smal a current, so that equation (23) is 
effectively the same as equation (22). A more direct connection appears 
if equation (22) is written as 


M' = 

E 0 a 



(22a) 


( See appendix C • ) 


MEASUREMENT PROCEDURES FOR VELOCITY 
AND TEMPERATURE FLUCTUATIONS 
Single Normal Wire 


It has teen seen that the. response of a hot-wire to simultaneous 
velocity and temperature fluctuations is given approximately ty 


e = 


®a " ®r 


IB. 


2 _ (R - R a ) 2 

1 R - aZn 2_ 

(R 0 a ) 2 


R - R„r 
— r- 5 i 

2iR„ L 


i 2 R - A(R - R a ) ~ 
J U 


(19) 


This can be written simply as 


e = 


4 + $ 


(19a) 


where 0 and 6 are functions of the mean-velocity calibration constants 
and of the hot-wire operating conditions. 


First consider the quantities that can be obtained from the root -mean- 
square output of single wires set perpendicular to the mean flow: 



( 24 ) 
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Since P and 5 do not increase in the- same proportion for a given 
increase in hot-wire temperature , readings of e^ at three different 
currents at the same point in the flow will give simultaneous linear 
algebraic equations from, which the three unknown quantities can he obtained: 



It immediately follows that 
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In general, because of the rather complicated forms of 3 and 6, 
these' determinants will be tedious to evaluate. However, the procedure 
may be made a bit more convenient by fixing the ratio RjjRg:!^ and setting 

up appropriate calculation charts. 


Possible Short Cuts 


In order to measure either 



to make either 6 or 3 negligible; that is. 


directly, it is necessary 




i^R - aZn 


(K - Ra) 

(Ro a ) 2 


2 


» 

or 

« 



The meaning of this can be found in terms of operating conditions at a 
given point in the flow by replacing the electrical quantities with fluid 
mechanical quantities as far as possible. Then the inequality can be 
written, after some algebra, as 




_ » 
0 or 
« 



( 26 ) 


For the entire range of practicable operation, the first term in 
the parenthesis in equation (26) Is larger than the second. Therefore, it 
is evident that,^for a given flow, 3 is maximized relative to 6 by 
minimizing R - R a , and the converse is also true. Whether either of 

these extremes can, in practice, make 3 » & or 3 « 8 can be checked 
numerically from a typical hot-wire calibration. It will be noted that 
the first term in the parenthesis is the only one of the three that can 
be varied by changing the hot-wire current. It turns out that, when 
both 6 and U are large, the only feasible inequality is that obtained 
by reducing the heating current so that 3 » 5; this ordinarily leads to 
the additional simplification, 

3 = i(Ra " 

for the sensitivity of the simple resistance thermometer. 

Equation (26) also gives the obvious result that reducing flT toward 
zero will eventually make 3 « 5, and then the wire is a simple anemometer.. 
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Another possibility for direct measurement is the use of two close 
parallel wires of different diameter, the voltage difference being the 
signal to the amplifier. Then, if 3-^ = 02 > the difference will 


measure 



if 


62, the difference will measure 



Of course, they must be far enough apart to be insensitive to lateral 
velocity fluctuations*. 


The general requirements for 3q = 02 a 1 *© extremely complicated if 

the general expression for 3 is used- However, for fairly rana l 1 currents, 
the second term in 3 is negligible compared with the first. In cases 
where this is true, 3]_ = 02 means 

h^(%-B ri )=l 2 ^(S a2 - Rr2 ) 


But _ _ 

®aq. ” ^r^ ^ a 2 


so that 

i^ ^i^ ( 27 ) 


This relation can be satisfied physically, but the problem of correct 
compensation for the two wires of different diameters would probably 
require that the fluctuating voltages be subtracted after leaving matched 
amplifiers, before entering the output circuit. 

The requirement for S 1 = & 2 is rather complicated, and this case 
will not ordinarily be necessary, since it is usually possible to make 

6 « 3- 


Directionally Sensitive Meters (e.g. , X and || ) 


In an ideal meter, the two wires are identical, so that the voltage 
fluctuations will be 0 

e, = <£+*£ 

1 0 u u 


d -3 c U V 

eo = 3~ + 5— - e— 
9 U U 


( 28 ) 
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Subtracting the voltages before amplification and measuring the mean 
square leads to 


v^2 (©i " © 2 )“ 

JJj 




(29) 


thus giving the lateral turbulence level. 

Taking the mean square output of each wire separately, there results 



After subtraction, 

4*11 (, 0 ) 

eu u 2 . 

The turbulent momentum and heat transfer are determined by repeating 
equation (30) at a different current and solving the resulting pair of 
simultaneous linear equations. Of course, the two terms in equation (30) 
are proportional to the turbulent shear and heat transfer only if the 
density fluctuations can be neglected (appendix D) . 

If the hot-wire instrument is unsymmetrical, the analysis is more 
complicated but essentially the same. For purely kinematic measurements, 
the generalization is given in reference 19* 

Just as in the case of the single wire, direct measurement of 

uv #U 

or -g is possible if 5« 0 or S» 3, respectively. Again it Is 

to be expected that in a fairly hot jet only the former condition is 
obtainable . 

In the general case, corresponding to equation (30), and in all other 
correlation measurements involving thermal and kinematic variables 
simultaneously , the direct measurement of correlation coefficients 
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^©•g* > 
direct 


^ 7 ' ^77 > etc *) is apparently 
determination is common practice 


not usually possible. Such 
in purely kinematic measurements. 


Correlations between Different Points 

For two single wires set parallel to each other and perpendicular 
to the. flow , assuming the ideal case of identical wires at the same 
sensitivity, the following equations apply: 


fl l u l 
e-. = 3 = + 

1 9 U 


t 3 2 hg 

3 2 = + 6r= 

d 9 U 


\ 


> 

J 


Adding and subtracting the fluctuating voltages gives 



(3D 


Expanding the foregoing expressions and subtracting gives 



and the three quantities in parentheses on the right side can be obtained 
by making readings at three different hot-wire currents. The physical 
significance of the cross-product term is not immediately evident. 

Of course, with 3 » 6, the thermal correlation function, and hence 
the thermal scale, is measurable directly. 
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The correlation between v fluctuations at different points can he 
determined by using two double wire meters, measuring the mean square of 
the sum of the differences of the individual meter differences. It can 
also be measured with sy mme tric ally inclined single wires. 


Spectrum 

'Q * ^ ' 

If ==■ and •=■ are the same order of magnitude in a flow with 

0 u 

appreciable velocity and temperature differences, it may be possible to 
measure only the spectrum directly. This is true for cases in which 

it is not feasible to make 6 because of wire temperature limitations. 

Then it may be necessary to obtain the u 1 spectrum with a meter made up 
of two wires of differing diameter as described previously. 

For temperature differences up to perhaps 100° C, however, it does 
not seem likely that the nature of the turbulence will be changed essentially 
in a flow of moderate velocity. Therefore, to a good approximation, it 
would seem satisfactory to determine the u 1 spectrum at reduced &■ 

It should be noted that the v 1 spectrum is directly measurable 
in all cases. 


MEASUREMENT OF MEAN VELOCITY AND CONCENTRATION 
IN AN ISOTHERMAL FLOWING MIXTURE OF TWO GASES 


As in the section dealing with the measurement of mean velocity and 
temperature, start with King's equation for thermal equilibrium. 


1 % 

R - R t 


= A + B 


ft 


where 


and 
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Again, A and B are not constants but depend upon the local physical 
constants of the flow and, therefore, are a function of the relative 
concentration. 

As indicated in the symbols, the concentration used in this analysis 
is the molal concentration. Further remarks on this are included in 
appendix E. 

The expressions for the mean physical constants of the mixture are 
written as 


P = P a + r Ap 


k = k a + f(r)Ak 


The thermal conductivity of gas mixtures is generally not a simple 
weighted average (reference 24) and varies considerably for different 
pairs of gases. If k(r) is a known function, the procedure is relatively 
simple, and, for most pairs of gases for which data are available (refer- 
ence 25 ) , f(r) can be well approximated by a quadratic or cubic poly- 
nomial. If k(r) is unknown, which is still the case for many pairs 
of common gases, it can be determined as one of the unknowns by an extra 
measurement at each point in the flow. Both cases are treated here. 



k(r) a Known Function 

For the discussion of k(r) as a known function, assume that f(r) 
can be approximated by a parabola: 


k = k + (|T.+ £r 2 )Ak 


(33a) 


Equation (2) can be written as 


_ i 2 R 
R - Re 


- igi[ k a + ( sr + 


bl 

R 


1/2 


(c Pa + r Ac^ l/2 (p a + r Ap) l/2 k a + (ir + cf 2 )^ u 1 / 2 ( 34 ) 
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Multiplying the first equation by N 2 and the second by N-j_ and 

subtracting, a solution can be obtained for |^k a + ^|F + ([r^Akj, and the 

resulting quadratic equation can be solved for F. Provided that k(T) 
is a monotonic function, only one of the two solutions of this quadratic 
will have any physical meaning. _This value of F can then be substituted 
into either equation to obtain U • This latter step can be done graphically 
if a mean-velocity calibration is drawn up of either hot wire; that 

•j I 

is, ^ — against \) U with straight lines of constant F. The two 

limiting lines r = 0 and f = 1.0 can be determined experimentally. 


For research in turbulent material diffusion it will be of particular 
interest to use a gas such that p g = p a (e.g., ethylene). In that case, 

a somewhat shorter technique can be used: A pitot tube will give tJ 

directly and then T can be determined from a single hot-wire reading 
and the family of calibration curves Just described. 


k(f) ‘an Unknown Function 


In the case of k(r) as an unknown function, it is convenient to 
work with the readings of two different hot-wires and a pitot tube at 
each point. 


*1 

(f>2 


The two hot-wires give 

“ M l[*a ♦ + Ni(c Pa ♦ f Ac p ) l/2 (p a ♦ r^[k a 

= H>[k a + f(DAk] + N 2 (c pa + F Ac p ) l/2 (p a ♦ f Ap) l/2 [k a 


1/2 ^ 
f(r)AkJ U'l/2 


-1 1/2 / 

f(r)AkJ u 1 / 2 


> (35 


where 




i g E 

E - E a 


M = 


aZ 

E 0 a 


N = 


blsfd 

R 0 a 
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Multiplying the first equation by N 2 and the second by N 
subtracting yields 

- 0U.. 

k a + f(r)Ak = — = K 

a MiN 2 - M^ a 


Substituting equation (36) into the expression for 


h = M i R a + N i( c P + r Ac p ) l/2 ( Pa + r Ap) l/2 Ka l/2 u l/2 


or 


(■ 


°Pa + r 


* v. 1/te 


A reading at the same point with a pitot tube gives 

q = ~(p a + f Ap)!! 2 


Thus, 


(p a + f Ap) 


1/2 (2 q ) 1 /2 


U 


which is substituted into equation (37) to give r(u): 


F ■ sbte® ' V 


However, equation (38) also gives f(u) directly: 


P _ l/2q 
- Apt ^2 p a 


and 

(36) 


(37) 


(38) 


(39) 


(38a) 
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And these two expressions are equated to give a cubic for U: 


_ 3 + 2q_A0p/Pa _ ^Pa\- 2 . 4^ = q 

K b 2 x^> Ac p/ K ^ 2 Ap 


m 


With the value of U at a point, T is obtained from equation ( 38 a) 
or equation (39)* The last unknown is the function k(f). 


The values of at several points immediately permit solution 

for f(T). For example, suppose it is assumed that f can be approximated 
— —2 -3 

by a cubic f = + agT + a^r . From readings at three different points 

in the flow, equation ( 36 ) yields three simultaneous equations for the 
unknown coefficients, 



+ r 2 a + r 3. 

+ 1 a 2 + 1 a 3 



_ 2 _ , ^ 

^ 2 a 1 + r 2 2 a 2 + r 3 


- k E 


2 a 3 - 


Ak 


> 


( 41 ) 



+ 



r 3 a 
3 3 




Of course, the three points should be rather widely spaced, and errors 
can be reduced by choosing several sets of three points at random. 



NACA TN No. 1864 


33 


RESPONSE OF A HOT-WIRE IN AN ISOTHERMAL FLOWING MIXTURE 
WITH VELOCITY AND CONCENTRATION FLUCTUATIONS 
Voltage Fluctuation 


Start again with King's equation 


i g R 
R - R c 


= — k + — (iCppku) 
R n a R n a\ P S 


, 1/2 


' 0 “ ^ 

and introduce small fluctuations in velocity and concentration: 


R = R + r 
k = k + k' 


C P = C P + C P 


p = p + p' 

U = U + u 

r = F + 7 


m 


From equation (33) > 


p' = /Ap 
c t/ = / Aci 




> 


( 43 ) 


Substituting equation ( 38 ) into King's equation, using appropriate 
binomial expansions, neglecting products of small quantities, and sub- 
tracting the averaged form of King's equation yields 


. (R - R a ) 2 




(44) 
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But, e =. ir; B can Be replaced according to the averaged equation; 
and p', Cp', and- 4' are given in equation (43). Thus, 


e 


(R - R a ) 2 — df Ak 1/ i 2 R 
lR a A "dr k + 2 \r - R a 


(R - E a ) 2 / i 2 ^ X 

2iR a - R a ‘ Jv 



(45) 


A more convenient form is obtained by substituting 



, ac p r 

c Pa + r ^ 


P 


1 + - 


'Pa 


T Ac, 


p 


l 


l + 


Pa 

f Ap 



Ak 

Also, the -=-terms can be combined. Then the final form of the equation 

giving the instantaneous voltage fluctuation due to concentration and 
velocity fluctuations is 



— y 

For r = 0 or «=0, equation (46) reduces to equation (3)> the 
response equation for the simple anemometer. In general, for concentration 
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and velocity fluctuations of the same order of magnitude, none of the 
terms in equation (46) is negligible. Typical values are given in 
appendix B. 

As mentioned in the previous section, investigations on material 
diffusion in turbulent flow will first be concerned with the mixing of 
two gases of equal densities, in which case Ap = 0. 


Time Constant 

From the section CONVENTIONAL HOT -WIRE -ANEMOMETER RESPONSE EQUATIONS, 

f - i 2 B - (A - B \TO)(B - E a ) (la) 

With equation (2), this again gives 

4.2ms dR i^R a (^e ” ^ 

R 0 a dt E e " R a 


which is equivalent to equation (4) of reference 2. Hence, with the usual 
restriction of small fluctuations, the same approximation is used in the 
denominator of the right side: 


R e ' R a = (R - R a ) + (B e " B) = R - R a 


so that equation (4) he comes approximately 


l+.2ms 

R 0 a 


dR 

dt 


i 2 R a (R e 

R - R 


- R) 
a 


or 


^(R - R) + 1(R - R e ) = 0 
dt m e 


M = 


4-2ms(R - R a ) 

i 2 R R n a 
a 0 


where 
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just as for the simple anemometer. Of course, the numerical value of the 
time constant in this case is a function of "both local mean velocity and 
mean concentration. 


MEASUREMENT PROCEDURES FOR VELOCITY AND 
CONCENTRATION FLUCTUATIONS 
Single Normal Wire 


The voltage fluctuation for a simple hot-wire set normal to the mean 
flow direction (equation (46)) may he written as 

6 = ^ f + ^ (^6a) 

The general procedure to he followed in using equation (46) to 

determine ffi ■ (r) , and (S) is similar to that for the temperature- 
fluctuation problem, with an essential difference; namely, in this case 

p 

the three readings of e taken at a given point in the flow cannot 
ordinarly he made simply with the same hot-wire operating at different 
temperatures. -This is due, of course, to the fact that p and 5 vary 
proportionally with changes in wire current (or temperature) at a given point 
in the flow. Instead, three wires with nonproportional mean calibration 
constants must he used. According to King's equation, this is possible 
only with three wires of differing diameter. However, in actual practice,: 
the ordinary differences in A and B between two wires of the same 
nominal diameter may provide sufficient variation for this method. 

It should also he pointed out that, since A is a function of the 
radiation loss, it may he possible to make at least two of the three 
necessary readings by using one hot-wire at widely differing temperatures. 

Thus, the equations to he solved are completely analogous to 
equation (25) and need not' he written out. 


Possible Short Cuts 

In general, it is not possible to make p » 6 or p « & for 
ordinary ranges of concentration and velocity. However, it may he possible 
to use two close parallel wires, the voltage of which is subtracted before 
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entering the amplifier- Then, if uq = \i 2 , the difference will 

/ \2 2 
measure , while, if 6q = 52, the difference will measure . 


As pointed out under MEASUREMENT PROCEDURES FOR VELOCITY AND TEMPERATURE 
FLUCTUATIONS, the wires must, of course, he far enough apart to have 
negligible sensitivity to v 1 - 


The conditions necessary for making uq = U 2 without simultaneously 
making 5j_ = 5 2 are so complex as to seem impracticable- The inverse 
condition is a bit simpler: 



(^ 7 ) 


But even this is extremely difficult to employ. In fact, perhaps the 
simplest conceivable assumptions which may still be satisfied by a pair 
of wires may be considered: Let B}^ = B 2 and R fll = R ag with the 

specific condition that A 1 4 ^2., since equal A's together with the 

other conditions would also make uq = Ug* Now equation (47) bec ome s 


h _ (*q - * a > 2 
i 2 (R 2 - R a ) 2 


(48) 


In order to predict the necessary currents, R^ and R2 must be 
eliminated- This leads to 



which should be solved explicitly for i2(iq). This is not generally 
possible. Therefore, if this approach is to be used for direct measure- 


ment of 



perhaps a feasible technique would be to satisfy equation (48) 
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by trial and error. Also, it should be noted that, since this assumes 
nearly identical wires operating at different temperatures, the time 
constants will be different, and separate matched amplifiers will be 
necessary. 


Directionally Sensitive Meters 

The general approach and procedure is completely analogous to that 
described for simultaneous temperature and velocity fluctuations. 

Assuming a symmetrical meter. 


7 _u v 

*1 ' T + Vs 


7 u v 

e 2 = H's + - er 

r u u 




> 




(49) 



Thus the turbulent shear and lateral diffusion of material are obtained 
by repeating equation (51) with an appropriately different instrument 
and solving the resulting pair of simultaneous equations. Of course, the 
terms in equation (51) are proportional to the turbulent shear and material 
transfer only if the density fluctuations can be neglected (appendix D) . 


Correlations between Different Points 

An approach identical with that for the temperature and velocity 
fluctuations gives, instead of equation ( 32 ), 
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and the three quantities in parentheses on the right side can he obtained 
by repeating the readings with two other appropriately different pairs of 
hot wires. 

Again, unlike the case in the section MEASUREMENT PROCEDURES FOR 
VELOCITY AND TEMPERATURE FLUCTUATIONS, neither nor can be 

r u 

measured directly. 

The correlation between v fluctuations at two different points 
is determinable as outlined In the aforementioned section. 


Spectrum. 

Since, in general, the value of p cannot be made either very large 
or very small compared to the value of 5, no direct measurement of either 
the y' spectrum or the u' spectrum can be made. 

In the case of gases whose densities and viscosities are not too 
widely different, the u' spectrum in the mixing gases is probably 
quite close to that in a similar flow of either component alone, for 
which case it is easily measurable* A comparison of this u 1 spectrum 
in the pure gas with an e' spectrum measured in the mixture may, at 
least, give a rough idea of the direction in which the y' spectrum 
varies as compared with u'. 

Again it should be noted that the v 1 spectrum is measurable 
directly in a flowing gas mixture with concentration fluctuations. 


CONCLUDING DISCUSSION 


In general, the errors Inherent in these new applications of the 
hot-wire anemometer are the same as for the simple anemometer (references 2 
26, 27 , and others) and hence need not be discussed here. 

Perhaps the only additional source of error arises simply from the 
added complication of many of these measurements. This Influences the - 
result in two ways: (a) The final numerical result Is obtained from 

more separate readings of fluctuating quantities than is the case in 
ordinary anemometric applications and (b) the appreciably greater time 
required for measurement at each point in the flow results in greater 
difficulty in maintaining the steady-state condition of the apparatus, 
both aerodynamic and electrical. 
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. It should also be noted that the hot-wire length corrections will 
be different for the different quantities, since each case must be 
computed from the correlation function of the physical quantity in 
question. 

The details of the present analysis assume the validity of King's 
equation in its conventional form. The accuracy with which this equation 
predicts the influence of velocity changes has been verified countless 
times; preliminary attempts to verify or disprove the rest of thi 3 
equation (appendix A) have been indecisive. Considerably more careful 
measurements must be made, both in different pure gases and in a pure 
gas at various temperatures. However, the general method presented here 
is equally applicable to a possibly more accurate heat -transfer equation 
that may be found in the future. 

The possible usefulness of a method for such detailed investigation 
of nonkinematic quantities in a turbulent flow is clear. For basic research, 
most of the measurable statistical quantities should be of interest; for 
some practical problems, perhaps the direct measurement of mean distri- 
butions and of turbulent transfer will be paramount. 

In principle, there appears to be no necessity for restriction of 
the method to measurement in a flow with only two simultaneous fluctuating 
quantities- However, the complication that arises upon consideration of 
simultaneous fluctuations in velocity, temperature, and concentration 
would appear to render this generalization impracticable. The idea of 
being able to make such measurements in a turbulent combustion zone is 
naturally of interest. However, this extremely important problem presents 
at least two certain major difficulties: (a) In reasonable density ranges, 

the flame temperatures are excessive and (b) the intermediate stages of 
most chemical reactions would involve dealing with- a mixture of several 
components instead of Just two. In addition to these difficulties, the 
possible catalytic action of the instrument in the reaction zone might 
also be noted as well as the probability that the hot-wire might act as 
an igniter. 

In casting about for possible methods of making direct measurement 
of the various kinematic quantities in a flowing gas mixture, the 
possibility should not be ignored of finding a pair of gases whose 
densities, specific heats, and thermal conductivities are so related that 
the net response to a concentration fluctuation is negligible. It can 
be seen from the numerical values in appendix B that such a situation is 
not beyond the realm of possibility. 


Calif omia . Institute of Technology 
Pasadena, Calif., June 10 , 19I+7 
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APPENDIX A 

REMARKS. ON KING ' S E QUATION 
General 


As pointed out by McAdams (reference 28) , the forced-convection term 
in King's equation can be written in terms of the product of a Reynolds 
number and a Prandtl number, so that the rate of heat loss from the wire 
is 


i 2 R = k(T - T a ) (a 1 + b' pc) (Al) 

where R under the radical is Reynolds number. 

■Inserting an extra d and introducing a heat- transfer coefficient h 

Vi rl 

and the Nusselt number T| = — McAdams writes this in the form 

k 7 

t) = *1 + KgNflc (A2) 

where R is Reynolds number. 

Conventionally, the first term in King's equation is described as 
including radiation and free convection, while the second term is due to 
forced convection. The latter appears to be quite reasonable, but it is 
apparent that free- convection effects, if appreciable, cannot enter the 
heat- transfer equation simply additivelyj the free-convection phenomenon 
is complex, and in most current applications the direction of buoyancy- 
induced velocity is perpendicular to the main flow. The form of A, as 
given by King, gives no clue as to its physical origin, although it is 
physically obvious that radiation must be included. The question is 
‘perhaps further clouded by the fact that measurements of the zero-velocity 
point in constant-resistance hot-wire calibrations seem to agree rather 
well with the A obtained by the extrapolation of the points measured at 
relatively large flow velocities. 

It would seem that, in view of the extended possibilities of the 
hot-wire anemometer as described in the present report, a complete 
reexamination of the fundamental heat- transfer equation for the flow 
past a heated cylinder at sma.11 Reynolds numbers is in order, with par- 
ticular emphasis on explicit inclusion of all the significant physical 
quantities. For example, if A actually does represent a free-convection 
term, it should include a Grashof number, which generally has a rather 
complicated variation with temperature. 

For the purposes of the present work, preliminary attempts have been 
made to check the validity of King's equation. These are described 
briefly in the following paragraphs . 
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Mean Calibration of Hot Wire in Air at Various Temperatures 

Because of the difficulty of obtaining a series of velocities 
at each of several temperatures with the present equipment, single 
measurements were made at various velocities and temperatures, 
keeping R - E a = Constant. Since King’s B is relatively insensitive 
to temperature changes and since this forced-convection term seems 
physically sound, the measured points were used to compute the inter- 
cept A on the assumption that the theoretical variation of B(T a ) is 
correct. Thus, an A was computed from each measured point instead 
of from a complete faired calibration curve at each temperature. In 
this computation, the data of Taylor and Johnston (reference 29) for the 
thermal conductivity of dry air were vised. The air was, in fact, fairly 
humid, but here only the relative values are important, and an error in 
absolute magnitude of k would have no effect. The specific heat was 
assumed constant (cp = 0.24) over the temperature range, in accordance 
with the International Critical Tables (reference 23) • For the density 
variation, the air was assumed to be a perfect gas. 

The results (fig- 3) show considerable scatter because this form 
of presentation involves small differences between two relatively large 
quantities, one of which has appreciable scatter itself. From the figure 
it appears that A(T a ) probably varies more slowly than k(T a )j however, 
the large amount of scatter precludes the deduction of an alternative, 
empirical expression. The agreement appears to be sufficiently good to 
„ justify the use of King's equation in the present, fairly rough measure- 
ments of velocity and temperature fluctuations in a free jet. 


Comparison of Mean Calibrations in Air and Carbon Dioxide 

In a further attempt to check King's equation, a hot-wire was 
calibrated in both air and carbon dioxide . Then, with the assumption 
that the forms of King's A and B are correct, the carbon dioxide 
calibration was computed from the air calibration. The curves are given 
in figure 4. 

o 

For the computation of Aco 2 and Bco 2 from A air and B air , the 

thermal conductivities, which correspond to those of pure gases, were 
taken from references 29 arid 30* The specific heats and densities were 
taken from reference 23, and both gases were assumed to be perfect gases. 

From figure 4 it can be seen that, although the agreement between 
measured and computed Bco 2 is good, the two values of Aqq^ are 

relatively far apart. 

The principal possible source of error in the determination of Aqq 2 
and BcOo from the air calibration is the use of thermal conductivities 
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measured in relatively, pure gases. The laboratory air undoubtedly con- 
tained appreciable amounts of carbon dioxide and water vapor (the measure- 
ments were made on a cloudy day); the jet carbon dioxide was a commercial 

type, stated as 99 | percent pure, probably containing air, water, and 

other impurities. Small amounts of impurities may have an appreciable 
effect upon thermal conductivity, much larger, for example, than the 
effect upon density and specific heat. Unfortunately, insufficient ■ 
experimental data exist on the thermal conductivity of gas mixtures, and 
the theoretical approach has proved to be tremendously complicated, even 
the first approximation for monatomic gases (reference 24). Thus, it 
would be extremely difficult to estimate the error due to impurities, 
even if a quantitative analysis were available. It is concluded that 
the measurements should be repeated with pure gases or that the actual 
thermal conductivities should be measured . 
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APPENDIX B 

TYPICAL NUMERICAL VALUES OF TERMS IN RESPONSE EQUATIONS 
Temperature and. Velocity Fluctuations 

The complete equation (l8) can he written in the form 



It is of interest to compare the relative magnitudes of the four terms 
as indicated. 

Assume the following typical numerical values: 

U = 1000 cm/sec d = 0.0006ll cm 

9 = 100° C Rq = 11.00 

T r = 290° abs . I = 0.30 cm 

i = 0.040 amp 

-S- = L = 0 .20 

^ ® a = 0.0037 (for platinum wire) 

A = 0.0035 
B = 0.00010 

k = 7 *55 x 10 ^ (reference 29) 
p = 7*00 X 10~ k 
= 0.25 

n = 1.62 x 10-7 (reference 29) 
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These values give 

E r = Eo[l + a(T r - 273)] = 11.69n 
Ea = E 0 [l + a(T a - 273)] = 15-7^ 
and, from the equation for thermal equilibrium, 

i*"E 

E = Ea + r p f= -- 20.74n 

A - + B VU 

The following numerical values are now obtained for the terms in 
equation (Bl) : 

(a) = 10.90 X lO"^ 

(b) = 1.22 X 10 -1 ^ 

(c) = -0.13 X 10" 1 *’ 

(d) = 3.16 X 10" 1 * 

Term (c) is only a little more than 1 percent of tern (a) and is 
apparently negligible for all anticipated operating conditions. 

Concentration and Velocity Fluctuations in 

Mixtures of Air and Carbon Dioxide 

Since there is no information available on precise measurements of 
the thermal conductivity of mixtures of air and carbon dioxide, it is 
assumed for the purpose s_of this calculation that the variation is 
linear) that is, f(T) = r. Then equation (46) can be written in the 
form 



(a) (b) 
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Assume the following typical numerical values: 

U = 1000 cm/sec Z = 0*30 cm p a = 1.165 x 10" 3 

r = °*5° i = 0.040 amp p = 1.785 x 10” 3 

O 

T r = T a = 290 abs • a = 0.0037 Cp a = 0.25 

7j" = ~y = ® ^ = 0.0020 Cp = 0.199 

8 

d = 0.0006ll cm B = 0.00010 Ak = -2.13 X 10"5 

3 0 = n - 0fi 1% = 5.93 X 10” 5 Ap = 0.620 x 10' 3 

kg = 3.8O X 10~5 ACp = -0.051 

These give 

(a) = -3.14 X 10” 4 

(b) = -0.72 X 10' 4 

(c) = 1.33 X 10" 4 

(d) = 3*16 x lo" 4 

and none of the terms is negligible . 
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APPENDIX C 

COMPENSATION RESISTANCE FOR SIMPLE RESISTANCE THERMOMETER 


With appreciable heating current, the ordinary anemometric expression 
for the time constant can he used: 


4 .2ms (R ~ R a ) 
E 0 a i 2 R Q 


which is more convenient than the equivalent fora. 


(Cl) 


M t = 4.2ms 1 

E 0 a (A - i 2 + B ft!) 

The compensation resistance is then given by 


(C2) 



(C 3) 


where L is the inductance of. the coil in the compensation network. 


However, in order to measure temperature fluctuation level directly 
it is convenient to operate with negligible heating current. Then (R - 
becomes too difficult to measure in a turbulent flow, and the limit must 
be reached with equation (C2) rather than with equation (Cl) . Then 


M ' = 1 ^ (cm 

( A + B \J U ) 

where, as in equation (C2), the values of A _and B to be used are the 
average values corresponding to temperature T a . 

In the actual process of measurement, perhaps the shortest method 

is to determine (a + B\[u)by a measurement of =3^- with the wire 

R - R a 

appreciably heated, since this reading is also necessary for the deterai- 
nation of mean velocity. 
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APPEND IX D 

TRANSFER IN TURBULENT SHEAR FLOW WITH DENSITY FLUCTUATIONS 


In the sections on measurement procedures, methods are given for 
the determination of shear, heat transfer, and material transfer by 
measurement with a directionally sensitive hot-wire meter. The procedures 
given actually lead to the quantities uv, "3v, and Tv. These, of course, 
are proportional to the respective transfer rates only in cases for which 
the density- and specific heat fluctuations are sufficiently small, in 
which cases the local mean density can he used to give 


Shear; 


T = -pUV 


Heat transfer; Q = -pCp-av 

Material transfer: D = -pyv 


(EL) 


Since the density fluctuations have been taken into account for the 
computation of the hot-wire response, it may be well to investigate the 
limitations introduced by the use of equation (Dl) . 


Temperature and Velocity Fluctuations 
The general expression for turbulent shear is 

e 

T = -^UV 

but 

— "'i 

P = P + p' 

U = U + u 

► 

V = v 
0 = 0 + ■d 


(T£) 


(D3) 


For turbulent shear flow with temperature gradient, it may reasonably be 

MV $ 

assumed that rr, ;=> and -=%■ are the same order of magnitude; in fact, in 


the previous discussion they are assumed small. It is of interest to 
find the restriction on 0 such that equation (Dl) may be applied with 
at least the same accuracy as may be expected from the measuring 
technique . 
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With equation (D3), equation (12) "becomes 

t = -puv - Up'v - p'uv 

and the third term is clearly negligible. 

Since the density fluctuation is due to temperature fluctuation, 
from equation (l7 ), 


P 


i 



Therefore, 


T 


puv 


+ £2 

T a 


0*) 


Hence, for equation (Dl) to he valid. 


uv » — dv 


or, 


uv dv e 
U 2 


(D5) 


If and are the same order of magnitude, this requires that 


U 2 


eu 


— « l (d6) 

T a 

If this condition is not satisfied with sufficient accuracy, the turbu- 
lent shear is calculable from equation (D4) after uv and dv have 
been obtained in accordance with equation (30) . 

The turbulent heat transfer is 

Q = -c pP 0Y (D7) 

but, with equation (D3) and with Cp = Constant, 


Q = -Cppdv - cp0p'v - Cpp'tfv 
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where the third term can he neglected. Substituting for p 1 yields 


__ g 

Q = -Cppdv + c p p =— fly 

•*a 


(D 8 ) 


so that, for equation (Dl) to he valid, it is required that 

~ « 1 (d6) 

T a 

>as before. If equation (d 6 ) is not satisfied, the turbulent heat 
transfer is easily calculable from 



after dv has been obtained from equation (30) . 


(D9) 


Concentration and Velocity Fluctuations 

With equation (D0), equation (D3), and the expression for the 
density fluctuation, 

p' = 7 Ap 

where T = T + 7 , there is obtained for the turbulent shear stress, 
approximately * 

T = -puv - U Ap 7 v (DIO) 


If equation (Dl) is to be reasonably accurate and if is roughly 

— . . V 2 

equal to the following inequality must be satisfied: 

ru 



(Dll) 


If equation (Dll) is not satisfied, uv and yv can be obtained from 
equation (51), and equation (DIO) gives the shearing stress o 

The turbulent material transfer is 


d = - P rv 


(D 12 ) 
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Introducing the mean values and fluctuations and neglecting the triple 
product of the small fluctuating quantities as before, equation (D12) 
becomes 


D = -p 7 v - rp'v 


(D13) 


But p' = y Ap, so that 

D = -(p + T Ap )yy (Dlh) 

and the condition for D in equation (Dl) to be valid is again 



(Dll) 


As before, equation (Dl4) may be used after ?v has been obtained from 
equation ( 51 ). 



52 


NACA TN No. 1864 


APPENDIX E 


MIXTURES OF PERFECT GASES 


In view of the sections discussing measurements in gaB mixtures, it 
seems appropriate to include a "brief review of pertinent relations 
occurring in the theory of mixtures of perfect gases. Extensive treat- 
ment is to "be found, for example, in reference 31 • 

The following notation is used: 

V total volume (for most problems in fluid mechanics, unit 

volume is considered) 

M-^ molecular weight of ith component 

Nj_ number of moles of ith component 

Pi partial pressure of ith component 

Pi density of ith component at pressure Pi 

Pi* density of ith component at pressure p 

R universal gas constant 

T absolute temperature 

p, p, M characteristics of mixture (M, apparent molecular weight) 

There follows immediately, 


P = YZ Pi 

i 


(El) 


N^ 
p i = — 


(E2) 


= Pi pj 


(E3) 
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Also, the equation of state is satisfied for each component and for 
the mixture : 


p ± V = N 1 ET 

(E4) 

pV = Y1 n ± et 

(E5) 

i 


p i 

(e6) 


and . 


P 



Since p = )> 

i 



so that 


(E7) 


M = 


i M i 


(E8) 


Eelative concentrations are of particular interest. Define the 
molal concentration of the ith component "by 



K 


(E9) 
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But, from, equations (E4) and (E 5 ) , 


Pi _ Ni 
P 

K 

so that 


pTI* 

II 

•H 

(E10) 

The following relation is of interest: 


p(Pi*, 7i) 


From equation (E3), 


.-C,. £(.,•*) 


or 



(Ell) 

For a mixture of two gases, 7g = 1 - y-j., and 


p = P2* + ^ Pl * - p 2 *) 7i 

(E12) 

A mass concentration can he defined as 


p i p i 

X± - r — = — 

YZ p k p 

K 

(E13) 

• 

A general expression for p(pp*, is rather 

equation (E13) can he transformed to 

complex. But 

X i = v~ “ 7i 

Z P K 

K 

(El4) 

Specializing at once to a two -component mixture, 
t Pi* 

i c = ~ y '- 

(E15) 
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"but, from equation (EI2), 


7l = 


P ~ P2* 


Pi* - P 2 * 

Substituting this into equation (E 15 ) and. solving for p yields 

P 1*P2* 


P = 


Pi* + ^P2* " Pl^^l 
which may be compared with equation (E12) . 


(El6) 


A relation between these two definitions of concentration is obtained 
by equating equations (E12) and (El6) : 


X.1 


Pl*?l 

P 2 * - ( Pl * - P 2 *)7 1 


(E17) 


or 


7i = 


P 2*H 

Pp* + ( p 2* - p l*)H 


(E18) 


The relative significance of these two concentrations in the field 
of fluid mechanics merits consideration. In investigations of mass 
transfer in gas flow, the mass concentration would appear to be the more 
significant parameter. For example, in a comparison of lateral rates of 
diffusion of momentum, heat, and material in a particular turbulent shear 
flow, the X, distribution is to be compared with the velocity and 
temperature distributions . 


On the other hand, a treatment of material diffusion prior to chemical 
reaction would obviously be more concerned with the molal concentration. 

For hot-wire measurements, the simplicity of equation (E12) relative 
to equation (El6) recommends the determination of molal concentration 
and fluctuations first, and, if necessary, the use of equation (E17)>to 
obtain mass concentration. ' 


Of course, if pp* = p2*, a case of particular interest in the 
investigation of material diffusion in turbulent flow, the distinction 
vanishes . 
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Figure 1.- Variation of hot-wire constants with air temperature. Physical 
constants from reference 23. R r , resistance at reference temperature: 
R a , resistance at local temperature. 









Figure 4.- Calibration of g- - mil platinum wire in carbon dioxide and air. Physical constants from 

references 23, 29, and 30. 




